Centrality dependence of elliptic flow and QGP viscosity 
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In the Israel-Stewart's theory of second order hydrodynamics, we have analysed the recent 
PHENIX data on charged particles elliptic flow in Au+Au collisions. PHENIX data demand more 
viscous fluid in peripheral collisions than in central collisions. Over a broad range of collision cen- 
trality (0-10%- 50-60%), viscosity to entropy ratio {rj/s) varies between 0-0.17. 
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I. INTRODUCTION 



One of the important results in Au+Au collisions at 
RHIC is the large elliptic flow in non-central collisions 
(H-Q. Large elliptic flows establish that in non-central 
Au+Au collisions, a collective QCD matter is created. 
Whether the matter can be characterized as the lat- 
tice QCD d i| predicted Quark-Gluon-Plasma (QGP) 
or not, is still a question of debate. Qualitatively, elliptic 
flow is naturally explained in a hydrodynamical model, 
rescattering of secondaries generates pressure and drives 
the subsequent collective motion. In non-central colli- 
sions, the reaction zone is asymmetric (almond shaped), 
pressure gradient is large in one direction and small in the 
other. The asymmetric pressure gradient generates the 
elliptic flow. As the fluid evolve and expands, asymme- 
try in the reaction zone decreases and a stage arise when 
reaction zone become symmetric and system no longer 
generate elliptic flow. Elliptic flow is early time phenom- 
ena. It is sensitive probe to, (i) degree of thermalisation, 
(ii) transport coefficient and (iii) equation of state of the 
early stage of the fluid 0-Q • 

Dissipative effects like viscosity reduce elliptic flow. 
The conversion of initial spatial anisotropy to momen- 
tum anisotropy is hindered in presence of viscosity and 
elliptic flow is reduced. QGP viscosity is an important 
parameter, theoretical estimates of the ratio, (shear) vis- 
cosity over the entropy density, rj/s cover a wide range, 

0. 0-1.0. String theory based models (ADS/CFT) give a 
lower bound on viscosity of any matter rj/s > 1/Att [loj . 
In a perturbative QCD, Arnold et al fill estimated rj/s ~ 

1. In a SU(3) gauge theory, Meyer [12j gave the upper 
bound rj/s <1.0, and his best estimate is i]/s=0. 134(33) 
at T = 1.165T C . At RHIC region, Nakamura and Sakai 
[13} estimated the viscosity of a hot gluon gas as ?7/s=0.1- 
0.4. Attempts have been made to estimate QGP viscosity 
directly from experimental data. Gavin and Abdel-Aziz 
[l4| proposed to measure viscosity from transverse mo- 
mentum fluctuations. From the existing data on Au+Au 
collisions, they estimated QGP viscosity as ?7/s=0.08- 
0.30. Experimental data on elliptic flow has also been 



used to estimate QGP viscosity. Elliptic flow scales with 
eccentricity. Departure from the scaling can be under- 
stood as due to off-equilibrium effect and utilised to esti- 
mate viscosity [l5j as, 7y/s=0.11-0.19. Experimental ob- 
servation that elliptic flow scales with transverse kinetic 
energy is also used to estimate QGP viscosity, rj/s ~ 0.09 
± 0.015 [lj], a value close to the ADS/CFT bound. From 
heavy quark energy loss, PHENIX collaboration [ljj es- 
timated QGP viscosity rj/s « 0.1-0.16. Recently, from 
analysis of RHIC data, in a viscous hydrodynamics , up- 
per bound to viscosity is given rj/s < 0.5 [H, [l9j]. Using 
a fully lattice based equation of state (lattice simulations 
for both the confined and the deconfined phase), we have 
estimated QGP visc osity as, rj/s = 0.15 ± 0.06 ^20]. In a 
more recent analysis [2l| , with a lattice based equation of 
state ( lattice simulations for the confined phase and non- 
interacting hadronic resonance gas for the hadron phase) , 
viscosity was estimated as, rj/s — 0.07 ± 0.03 ± 0.14. In 



21], estimate of viscosity was obtained from analysing 
STAR data on <f> mesons multiplicity, mean px and in- 
tegrated «2- In the same model from a study of 'scaling 
departure' of elliptic flow, viscosity over the entropy ra- 
tio was estimated as, rj/s — 0.12 ± 0.03 [2^]. In the 
present paper, we have studied the centrality depen- 
dence of elliptic flow in ideal and viscous hydrodynam- 
ics. From a direct fit to the PHENIX data [23[ on the 
charged particles elliptic flow in 0-60% centrality colli- 
sions, we have also obtained an estimate of QGP viscos- 
ity, rj/s = 0.115 + 0.005, very close to the value obtained 
from the study of scaling departure of elliptic flow [22j . 
Our analysis also indicate that the PHENIX data on the 
centrality dependence of differential elliptic flow require 
stronger viscosity in peripheral collisions than in central 
collisions. Over the broad range of collision centrality 
(0-10% - 50-60%), viscosity over the entropy ratio varies 
between -0.17. 
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The paper is organised as follows: in section II, we 
briefly explain the 2nd order Israel-Stewarts theory of 
dissipative hydrodynamics. Equation of state, initialisa- 
tion of the fluid is described in section III. Results are 
discussed in section VI. The summary and conclusions 
are given in section V. 
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II. HYDRODYNAMIC MODEL 

In the following, we consider a baryon free fluid with 
only shear viscosity. Bulk viscosity is neglected. In the 
Israel-Stewart's theory of 2nd order dissipative hydrody- 
namics, equation of motion of the fluid is obtained by 
solving, 



0, 



1 



Dtt^ = (tt"" - 2^V < 'V /> ) 

- [u"7r" A + u lJ ^ x ]Du x . 



(1) 



(2) 



Eq[T] is the conservation equation for the energy- 
momentum tensor, T^ v = (e + p)u^u v — pg^ u + n^ v , 
e, p and u being the energy density, pressure and fluid 
velocity respectively, is the shear stress tensor (we 
have neglected bulk viscosity and heat conduction) . Eqj2] 
is the relaxation equation for the shear stress tensor tt^ v . 
In Eql2l D — u^d^ is the convective time derivative, 
v <a^> = + V 1 ^) - \{d.u){g^ - u»u v ) is a 

symmetric traceless tensor. r\ is the shear viscosity and 
7v is the relaxation time. It may be mentioned that in 
a conformally symmetric fluid relaxation equation can 
contain additional terms (24j . 

Assuming boost-invariance, EqsJT] and [2] are solved in 
(r = y/t 2 — z 2 ,x,y,T) s = h In jz3) coordinates, with the 
code " AZHYDRO-KOLKATA" Z ', developed at the Cy- 
clotron Centre, Kolkata. Details of the code can be found 
in [Hi. Within 10% or less, AZHYDRO-KOLKATA sim- 
ulation reproduces Song and Heinz's [lij result for tem- 
poral evolution of momentum anisotropy e p . 
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FIG. 1: Lattice simulation for entropy density is compared 
with the same in model EOS, lattice+HRG. The filled circles 
are lattice simulation 6] for s/T z . The solid line is the same 
in model EOS (see text). 



simulation than in a HRG model. It is difficult to resolve 
whether the discrepancy is due to failure of HRG model 
at lower temperature or due to the difficulty in resolving 
low energy hadron spectrum on rather coarse lattice [f|. 

From the entropy density, using the thermodynamic 
relations, 



P(T) 
s(T) 



s(T')dT' 
Ts-p, 



(4) 
(5) 



III. EQUATION OF STATE 



pressure and energy density can be obtained. 



Eqs lll2l are closed with an equation of state (EOS) 
p = p{e). Lattice simulations [j| Q indicate that the 
confinement-deconfinement transition is a cross over, 
rather than a 1st or 2nd order phase transition. In FigfTJ 
a recent lattice simulation [6J for the entropy density is 
shown. We complement the lattice simulated EOS [a] by 
a hadronic resonance gas (HRG) EOS comprising all the 
resonances below mass 2.5 GeV. In Figfl] the solid line 
is the entropy density of the '"lattice +HRG" ' EOS. The 
entropy density is obtained as, 



IV. INITIALISATION OF THE FLUID 



0.5[1 — tanh(x)]s 



HRG 



0.5[1 + tanh(x)}siatuce (3) 



where siattice and shrg are entropy density from lattice 
simulations and HRG model, x = T ^ c ■ In the present 
simulation, we have used cross over temperature, T c =196 
MeV and AT = 0.1T C . Compared to lattice simulation, 
entropy density in HRG drops slowly at low tempera- 
ture. It is consistent with observation in [|| , that at low 
temperature, trace anamoly, drops faster in lattice 



Solution of partial differential equations (Eqs ll|2|) re- 
quires initial conditions, e.g. transverse profile of the 
energy density (s(x,y)), fluid velocity (v x (x,y),v y (x,y)) 
and shear stress tensor (^"(x, y)) at the initial time Tj. 
One also need to specify the viscosity (rj) and the relax- 
ation time (tv). A freeze-out prescription is also needed 
to convert the information about fluid energy density 
and velocity to particle spectra and compare with ex- 
periment. 

We assumed that the fluid is thermalised at Ti=0.6 
fm [26[ and the initial fluid velocity is zero, v x (x,y) = 
v y (x,y) = 0. Initial energy density is assumed to be 
distributed as 126] 



e(b, x, y) = e l [0.75N part (h, x, y) + 0.25N coU (b, x, y)], 

(6) 

where b is the impact parameter of the collision. N part 
and N co u are the transverse profile of the average num- 
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FIG. 2: (color online) Viscous suppression of elliptic flow in 
Au+Au collisions. The solid, dashed and mid-dashed lines are 
fractional suppression of elliptic flow in viscous fluid evolution 
with 77/s=0.08, 0.12 and 0.16 respectively. 



FIG. 3: (color online) The symbols are PHOBOS data on 
the centrality dependence of charged particles integrated V2- 
The solid, dashed, medium-dashed and short-dashed lines are 
simulated flow with rj/s=0, 0.08, 0.12 and 0.16 respectively. 



ber of participants and average number collisions re- 
spectively, calculated in a Glauber model. Central en- 
ergy density, £$ is a parameter and does not depend 
on the impact parameter of the collision. The shear 
stress tensor was initialised with boost-invariant value, 
tt xx — ir vv — 2rj/3Ti, ir xy =0. For the relaxation time, we 
used the Boltzmann estimate = 3ri/2p. Finally, the 
freeze-out was fixed at Tp=150 MeV [27j . 

TABLE I: Initial central energy density (si) and temperature 
(Ti) of the fluid in b=0 Au+Au collisions, for different values 
of viscosity to entropy ratio {ij/s). The predicted <f> meson 
multiplicity and mean pr are also noted. 



T]/s 


Si 

(GeV/fm 3 ) 


T, 
(MeV) 


dJV* 
dy 


< Pt > 
(GeV) 





35.5 ±5.0 


377.0 ± 13.7 


7.96 


1.02 


0.08 


29.1 ± 3.6 


359.1 ± 11.5 


8.01 


1.06 


0.12 


25.6 ±4.0 


348.0 ± 14.3 


8.22 


1.11 


0.16 


20.8 ± 2.7 


330.5 ±11.3 


8.13 


1.17 



Assuming that throughout the evolution, viscosity to 
entropy ratio (77/s) remains a constant, we have simu- 
lated Au±Au collisions for four values of rj/s, (i) rj/s=0 
(ideal fluid), (ii) 77/s = 1/4tt « 0.08 (ADS/CFT lower 
limit of viscosity) , (hi) ?y/s=0.12 and (iv) ?7/s=0.16. Note 
that during the evolution, the fluid cross over from QGP 
phase to hadronic phase. Constant 77/s can be thought 
over as a space-time averaged 77/s. As the entropy den- 
sity of the QGP phase is more than that of a HRG (see 
FigHJ), constancy of r/fs during the evolution implicitly 
assume t/qqp > t/hbg- It also fixes the temperature de- 
pendence of r], Variation of rj with temperature is the 
same as that of the entropy density. 



Only parameter left to be initialised is the central en- 
ergy density . We fix £j such that entropy at the freeze- 
out is the same either in ideal or in viscous evolution. 
This is done by reproducing the STAR data [28| on (j) 
meson multiplicity in 0-5% Au±Au collisions. In the 
present work, we have neglected resonance contribution. 
We choose <fi meson multiplicity as they are not affected 
by resonance decays. In table]!] initial central energy 
density (ei) and temperature (X^) required to fit STAR 
data on 4> meson multiplicity in 0-5% Au±Au collisions 
are noted. The error in Si or in Ti corresponds to statis- 
tical and systematic uncertainty in STAR measurements 
[28j ] . In viscous fluid evolution, entropy is generated. 
More viscous is the fluid, more is the entropy genera- 
tion. As a consequence, viscous fluid requires less ini- 
tial energy density (or temperature) than an ideal fluid. 
For example, compared to ideal fluid, in minimally vis- 
cous (77/5=0. 08) fluid, initial energy density is reduced 
by 18%. In fluid with viscosity 77/5=0. 16, reduction is 
even more, ^40%. The predicted central values of <fi me- 
son multiplicities and mean pj* are also shown in tabled] 
They should be compared with STAR measurements [281 ] , 

^ X = 7.95 ±0.74, and < p% >= 0.977 ±0.064 (statis- 
tical and systematic error included). Experimental data 
on (f> meson multiplicity and mean px in 0-5% central- 
ity Au±Au collisions are simultaneously explained only 
for rj/s <0.12. We have not shown it here, but hydro- 
dynamic evolution of the fluid initialised as stated here, 
reproduce experimental pr spectra of pions and kaons at 
Pt >1 GeV (pt < 1 GeV part of the spectra is under 
predicted due to neglect of resonance contribution). It 
also reproduces the px spectra of <fi mesons. The proton 
spectra however is underestimated by a factor of ~ 2. 




FIG. 4: (color online) Fractional suppression of differential V2 
in Au+Au collisions in the 0-60% centrality ranges of colli- 
sions. 



V. ELLIPTIC FLOW IN IDEAL AND VISCOUS 
HYDRODYNAMICS 

A. Centrality dependence of integrated V2 

Let us first study the effect of viscosity on pt in- 
tegrated elliptic flow. In Fig|2j fractional decrease in 
charged particles (pions, kaons and protons) elliptic flow, 
(1 — v% ls /v 2 d ), as a function of collision centrality is stud- 
ied. As mentioned earlier, we have neglected resonance 
decays. Resonance contribution reduces elliptic flow, 
mostly for low mass particles at low p T [2{|, e.g. in the 
Pt range, < pr < 1 GeV , V2 for pions is reduced by ~ 
0-30%. At pr > 1 GeV, resonance contribution to v<x is 
negligible [29} . Neglect of resonance contribution will in- 
crease integrated V2, however, the effect will be minimised 
in the ratio v 2 ls /v 2 d . In collisions with N part > 100, 
integrated vi is suppressed by ~5%, 10% and 15% for 
viscosity ?7/s=0.08, 0.12 and 0.16 respectively. There is 
also an indication that in peripheral collisions, V2 is more 
suppressed than in central collisions. In very peripheral 
collisions, N part <100, V2 in viscous fluid is more than 
that in ideal fluid. However, in very peripheral collisions, 
applicability of hydrodynamics is questionable [26| . We 
donot think that present simulations are reliable for pe- 
ripheral (N par t <100) collisions. 

In FigUH we have compared simulated (integrated) el- 



FIG. 5: (color online) In six panels, PHENIX measurements 
for elliptic flow in 0-10%, 10-20%, 20-30%, 30-40%, 40-50% 
and 50-60% Au+Au collisions are shown. The solid, dashed, 
medium dashed and short dashed lines are elliptic flow in 
hydrodynamic simulations with r]/s=0, 0.08, 0.12 and 0.16 
respectively. 

liptic flow with PHOBOS measurements [3(| for charged 
particles elliptic flow. In FigEl the solid, dashed, 
medium-dashed and short-dashed lines are the charged 
particles elliptic flow in the present simulation, with 
r]/s=0, 0.08, 0.12 and 0.16 respectively. Evidently, PHO- 
BOS data on the centrality dependence of i>2 prefer vis- 
cous fluid rather than ideal fluid. However, one also note 
that in central collisions N part > 250, simulated v 2 in 
ideal fluid evolution are more close to the experimental 
data than in viscous evolution. For N par t > 250, data 
are better explained in ideal dynamics than in viscous 
dynamics. Data at N part < 250 on the other hand are 
better explained in viscous evolution than in ideal evolu- 
tion. 



B. Centrality dependence of differential v-2 

Viscous suppression of differential V2 is studied in 
FiglU In FiglU in three panels, for ?y/s=0.08, 0.12 and 
0.16, fractional decrease in v 2 in 0-10%, 10-20%, 20-30%, 
30-40% and 40-50% Au+Au collisions are shown. In all 
the centrality ranges of collisions, viscous suppression is 
more at high pt than at low pt- Suppression is also cen- 
trality dependent, more in peripheral collisions than in 
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FIG. 6: x 2 / N for the PHENIX measurements of v 2 {2} in 
0-60% Au+Au collisions, as a function of rj/s are shown. The 
solid line is a fit to the \ 2 /N values by a parabola. 



central collisions. For example, for 77/ s= 0.12, at a fixed 
pr ~1 GeV, compared to a central collisions, in periph- 
eral collisions, v 2 suppression is increased by a factor of 
~ 2. At large pt however, viscous suppression tends to 
saturate for all collision centrality. 

In Figl5l we have compared the presently simulated v 2 
with PHENIX measurements [23 1. I n FigJS] colored sym- 
bols are PHENIX measurements 23| for charged particles 
elliptic flow in 0-10%, 10-20%, 20-30%, 30-40%, 40-50% 
and 50-60% Au+Au collisions. PHENIX collaboration 
measured charged particles v 2 upto pr ~ 8 GeV. In FigJSJ 
measurements upto pr=3 GeV are shown only. Hydro- 
dynamic models are not well suited for high px particles. 
In order to study non-flow effects that are not correlated 
with the reaction plane, as well as fluctuations of tj 2 , 
PHENIX collaboration obtained v 2 from two indepen- 
dent analysis, (i) event plane method from two indepen- 
dent subdetectors, v 2 {BBC} and v 2 {ZDC- SMD} and 
(ii) two particle cummulant i>2{2}. «2{2} from two parti- 
cle cummulant and V2{BBC} or v 2 {ZDC — BBC} from 
event plane methods agree within the systematic error. 
It may also be mentioned here that t>2{2} in PHENIX is 
lower than that 1*2 {2} in STAR measurements, but they 
agree within the systematic error. All the three measure- 
ments of v 2 are shown in FigJS] 

In Fig[5j the solid, dashed, medium dashed and short 
dashed lines are simulated elliptic flow in fluid evolution 
with (i) r]/s=0 (ideal fluid), (ii) 77/3=0.08, (hi) ?//s=0.12, 
and (iv) ^/s=0.16, respectively. Comparison of simu- 
lated elliptic flow with PHENIX measurements indicate 
that in central (0-10%) collisions, hydrodynamic evolu- 
tion produces less v 2 than in experiment. For example, 
at pt ~1 GeV, ideal fluid evolution underestimate of 
PHENIX measurement of v 2 {ZDC - SAID} by -40%. 
Viscosity reduces v 2 and in viscous evolution v 2 is even 
more under predicted. That in central collisions, ideal 



hydrodynamic model, with Glauber model initial condi- 
tion produces less v 2 than in experiment is well known 
(e.g. see Fig. 9 of Ref.[8]). In a hydrodynamic model, el- 
liptic flow depends on the initial spatial eccentricity. In a 
Glauber model, initial eccentricity in central collisions is 
small and elliptic flow is not developed. Initial spatial ec- 
centricity is comparatively large in Color Glass Conden- 
sate (CGC) model initial conditions [3l|. Hydrodynamic 
with CGC model initial condition correctly predict ellip- 
tic flow in central collisions [3l[ . Apparently CGC initial 
condition is better suited for central Au+Au collisions. 
However, since a hydrodynamic model have quite a large 
number of parameters, e.g. initial time, initial energy 
density, equation of state, freeze-out temperature etc, un- 
less a detailed study is made, it is difficult to conclude 
that hydrodynamic models with Glauber model initial 
condition do not correctly predict v 2 in central collisions. 

In mid-central or in peripheral collisions, elliptic flow 
is over predicted in ideal fluid evolution. Data are better 
explained in viscous fluid evolution. To obtain a quanti- 
tative idea about the fit obtained to the PHENIX data by 
the hydrodynamic simulations, we have computed \ 2 /N , 



X 2 /N 



1 

N 



i=N 

E 

»=i 



(EX(i)-TH(i))< 
ERR(i) 2 



(7) 



where EX{i) and ERR(i) are the PHENIX measure- 
ments for W2{2} and its error, TH(i) is the hydrody- 
namic simulations for v 2 . As mentioned earlier, we have 
neglected resonance production. Resonances contribute 
mainly at low pr, Pt < 1 GeV [29| . To remove the un- 
certainty due to neglect of resonance contribution, in the 
X 2 /N computation, we have included data only the in 
the p T range 1 GeV < p T < 3 GeV. In tableHH the 
X 2 /N values for different centrality ranges of collisions 
are noted. One note that in most central (0-10%) colli- 
sions, the minimum value is obtained in ideal fluid evolu- 
tion. More peripheral collisions prefer viscous fluid. It is 
also be noted that the minimum x 2 /A r in 0-10% central- 
ity collision is factor of 2 or more larger than minimum 
X 2 /N obtained in other collision centrality. Compared to 
peripheral collisions, elliptic flow in very central collisions 
are not well reproduced in a hydrodynamic model. 

X 2 /N for the combined data sets as a function of 77/s 
is shown in FigEJ The solid line in Figgis a parabolic 
fit to x 2 /N. Best fit to the PHENIX combined data sets 
is obtained for 77/s = 0.115 ± 0.005. The estimate is 
well within the upper bound of 77/s obtained in (l8l. Il9j|. 
The result is very close to the estimate rj/s = 0.12 ± 
0.03, also obtained in the same model, from the study of 
scaling violation of elliptic flow [22j|. It also agree with 
the estimate r?/s=0.11-0.19 0, 77/s = 0.09 + 0.015 
also obtained from the analysis of elliptic flow data. 

PHENIX measurements of charged particle elliptic 
flow in 0-60% Au+Au collisions, are best explained with 
viscosity to entropy ratio 77/s «0.115. However, as shown 
in Fig[Bl the minimum (x 2 /-^)mm ~18, is comparatively 
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TABLE II: x*/N for PHENIX v 2 {2} in 0-60% Au+Au colli- 
sions, in fluid evolution with viscosity to entropy ratio, r//s=0, 
0.08,0.12 and 0.16. 





X 2 /N 


coll. centrality 


77/s = 0.0 


r//s=0.08 


?j/s=0.12 


?j/s=0.16 


0-10% 


11.79 


17.87 


21.51 


25.27 


10-20% 


19.48 


5.78 


33.61 


93.77 


20-30% 


215.21 


5.38 


26.89 


152.98 


30-40% 


506.16 


51.53 


6.33 


121.56 


40-50% 


669.40 


125.31 


8.13 


37.70 


50-60% 


494.09 


140.72 


38.52 


2.75 



large (a good fit to data demand (% 2 /N) min ~1 ). Com- 
paratively large value (\ 2 /N) min «18, indicate that hy- 
drodynamic model with a fixed viscosity to entropy ra- 
tio do not predict 'accurately' centrality dependence of 
elliptic flow. Indeed, one observe from table|TT] that in 
more peripheral collisions, PHENIX data demand more 
viscous fluid. For example, in 0-10% Au+Au collisions, 
minimum y 2 /N is obtain in ideal fluid (r]/s=0) evolution. 
In 10-20% and 20-30% collisions, \ 2 / N is minimum for 
?7/s=0.08. In 30-40%, 40-50% centrality collisions, mini- 
mum X 1 /N is obtained for 77/s=0.12. In 50-60% central- 
ity collisions, minimum x 2 /iV is obtained at still higher 
value, ?7/s=0.16. Apparently, present analysis suggests 
that more viscous fluid is produced in peripheral Au+Au 
collisions than in central collisions. To obtain the central- 
ity dependence of 77/s, in each collision centrality, x 2 /N 
as a function of 77/s are fitted by a parabola. From the 
minimum of the parabola, best fitted 77/s is obtained. In 
table llll) best fitted 77/s as a function of collision cen- 
trality are noted. 77/s smoothly increases as the collision 
become more and more peripheral. While in central col- 
lisions, a nearly perfect fluid is produced, more viscous 
fluid is produced in peripheral collisions. 

Can we understand the centrality dependence of 77/s? 
Present paradigm is that 77/s has a minimum, possibly 
with a cusp, around the critical temperature T = T c 
[32j]. The centrality dependence of 77/s, as obtained in 
the present analysis, is not at variance with the prevail- 
ing paradigm. Rather it indicates the increasingly im- 
portant role of hardronic matter in the development of 
elliptic flow in peripheral collisions. As mentioned earlier, 
viscosity to entropy ratio, as obtained here, is averaged 
over the space-time. Both the QGP and the hadronic 
phase contribute to the average. Schematically, one can 
write, 

?7 77 ISP 77 had 

- = (1 - /had ) ( — ) (Tqgp) + /had(-) {Trad) 
s s s 

(8) 

where /had is the fraction of the hadronic matter, 
(f) 9ffP (Tqgp) is the viscosity of QGP matter at average 
temperature Tqgp and (^) haa (Th ad) is the viscosity of 
the hadronic matter at average temperature Thad- In 



TABLE III: Best fitted 77/s as a function of collision central- 
ity. Also shown are the fraction of hadronic matter (/had), 
spatially averaged temperature of the QGP matter ({T)qgp), 
and the hadronic matter ((T)had) at the initial time n. 



coll. centrality 


77/s 


J HAD 


(T)qgp 
(MeV) 


(T)had 
(MeV) 


0-10% 


+ 0.03 


0.20 


298.7 


173.0 


10-20% 


0.051 + 0.008 


0.27 


283.2 


172.8 


20-30% 


0.087 ± 0.004 


0.31 


267.9 


172.6 


30-40% 


0.109 + 0.003 


0.35 


254.7 


172.4 


40-50% 


0.134 + 0.004 


0.41 


239.6 


172.4 


50-60% 


0.169 + 0.005 


0.49 


222.7 


171.8 



table Hill we have noted /had at the initial time. In 
0-10% collisions, only ~ 20% of the matter is in the 
hadronic phase, the rest (~80%) is in the deconfmed (or 
QGP) phase. Fraction of the hadronic matter increases 
in peripheral collisions and in 50-60% centrality Au+Au 
collisions, at the initial time, ~ 50% of the total matter 
is hadronic. Even if elliptic flow develops early in the 
evolution, viscosity of hadronic matter will play impor- 
tant role in the development of the elliptic flow, more so 
in peripheral collisions. In table lllll we have also noted 
the spatially averaged initial temperature of the QGP 
phase (Tqgp) and the hadronic phase, (Thad)- As ex- 
pected, Tqgp decreases in peripheral collisions, spatially 
averaged initial temperature of the hadronic phase, on 
the other hand is nearly constant, independent of colli- 
sions centrality (the temperature range of the hadronic 
matter is limited, 150MeV < T HAD < 196MeVL Since 
(i]/s) c,9P decreases with decreasing temperature [32|, con- 
tribution of the QGP phase to the space-time averaged 
77/s will be less in a peripheral collision than in a central 
collisions. Since hadronic fraction increases in periph- 
eral collisions, contribution of the hadronic phase to the 
space-time averaged 77/s will increase in peripheral colli- 
sions. Centrality dependence of extracted 77/s can qual- 
itatively understood as due to increased contribution of 
hadronic phase and decreased contribution of the QGP 
phase in peripheral collisions than in a central collisions. 

The present estimate ?7/s=0-0.17 in 0-60% centrality 
range, must be treated with caution. In the present sim- 
ulations, we have neglected bulk viscosity. Experimental 
data, include the effect of bulk viscosity, if there is any. 
Neglect of bulk viscosity, will artificially increase the ef- 
fect of (shear) viscosity. In general, bulk viscosity is an 
order of magnitude smaller than shear viscosity. But in 
QGP, it is possible that near the cross-over temperature, 
bulk viscosity is large [33|, [34| . Effect of bulk viscosity 
on particle spectra and elliptic flow is studied in (35j |. It 
appears that even if small, bulk viscosity can have visible 
effect on particle spectra and elliptic flow. The present 
estimate then must be considered as an upper bound on 
QGP viscosity. Also, we have not considered any sys- 
tematic error in evaluation of 77/s due to our uncertain 



7 



knowledge about various parameters of the hydrodynam- 
ics model. Systematic error in hydrodynamic evaluation 
of viscosity could be large. Indeed, in [2l|, from a si- 
multaneous fit to <p meson multiplicity, mean px and in- 
tegrated «2, viscosity to entropy ratio was estimated as, 
Tj/s = 0.07 ±0.03 ±0.14, where the 1st error is statistical 
and the 2nd one is the systematic error. Large (~ 100%) 
systematic error arises due to uncertainty in initial time, 
initial energy density profile, initial fluid velocity, freeze- 
out condition, finite accuracy of computer code etc. Even 
then as noted in [2l[ the source of systematic error is not 
exhaustive. We expect the systematic error in the present 
evaluation of viscosity to entropy ratio will be of the same 
order ~ 100%. 



VI. SUMMARY AND CONCLUSIONS 

To summarise, we have studied effect of (shear) viscos- 
ity on elliptic flow. To obtain a meaningful comparison 
between flows in ideal and viscous dynamics, the fluid 
was initialised to reproduce (f> meson multiplicity in 0-5% 
Au±Au collisions. The initialisation ensures that irre- 
spective of fluid viscosity, the final state entropy is the 



same. Elliptic flow is suppressed in viscous fluid evolu- 
tion, more viscous is the fluid, more is the suppression. 
Depending on viscosity (?7/s=0.08-0.16), in central and 
mid-central (N part > 100 ) collisions, integrated i>2 is 
suppressed by 5-15% only. Suppression of differential v-i 
is pr dependent. t>2 (pt) is more suppressed at large pr 
than at low pr- Suppression of Vaipr) is also central- 
ity dependent, suppression is more in peripheral than in 
central collisions. Centrality dependence of elliptic flow 
suppression however reduces at large pt- 

We have also compared simulated flow in ideal and 
viscous dynamics with experiments. PHENIX data [23| 
on the differential elliptic flow in the centrality range 0- 
60% are best described in viscous fluid evolution with 
Tj/s = 0.12 ± 0.005. However, it was also indicated that 
the data demand more viscous fluid in more peripheral 
collisions. For example, in central collisions (0-20%), 
PHENIX data are best explained with small viscosity, 
rj/s « 0-0.05. In more peripheral collisions, e.g. in 40- 
60% collisions, data demand more viscous fluid, rj/s rs 
0.13-0.17. Centrality dependence of QGP viscosity can 
be understood as due to increased contribution of the 
hadronic matter in the development of elliptic flow in 
peripheral collisions. 
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